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Abstract. A method of finding relative periodic orbits for differential equations with continuous 
symmetries is described and its utility demonstrated by computing relative periodic solutions for the 
one-dimensional complex Ginzburg-Landau equation (CGLE) with periodic boundary conditions. A 
relative periodic solution is a solution that is periodic in time, up to a transformation by an element 
of the equation's symmetry group. With the method used, relative periodic solutions are represented 
by a space-time Fourier series modified to include the symmetry group element and are sought as 
solutions to a system of nonlinear algebraic equations for the Fourier coefficients, group element, 
and time period. The 77 relative periodic solutions found for the CGLE exhibit a wide variety of 
temporal dynamics, with the sum of their positive Lyapunov exponents varying from 5.19 to 60.35 
and their unstable dimensions from 3 to 8. Preliminary work indicates that weighted averages over 
the collection of relative periodic solutions accurately approximate the value of several functionals 
on typical trajectories. 
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1. Introduction. The complex Ginzburg-Landau equation (CGLE) is a widely 
studied partial differential equation with applications in many areas of science. (See 
44 and references therein.) It has become a model problem for the study of nonlinear 
evolution equations with chaotic spatio-temporal dynamics. In this paper we work 
with the CGLE with cubic nonlinearity in one spatial dimension, 



and, as a convention, when we refer to the CGLE we mean equation i|l.l|) with 
the boundary conditions l|1.2fl . unless otherwise noted. Equation (|1.1|) describes the 
evolution of a complex- valued field A(x,t). The parameters R, v, and fj, are real. 
When R < all solutions converge to A = in forward time, but when R > 
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(1.1) 



with periodic boundary conditions 



A(x,t) = A(x + 2n,t), 
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there is, in general, nontrivial spatio-temporal behavior. The parameters v and y, 
are measures of the linear and nonlinear dissipation, respectively. (See, for example, 
P ESI EH for further details.) It is known that in the form (JTTJl the CGLE 
generates a continuous semiflow on a variety of spaces |5*l l21ll3*9"ll56| . 

The CGLE has a three-parameter group T 2 x R of continuous symmetries gen- 
erated by space-time translations and a rotation of the complex field A. That is, if 
A(x, t) is a solution of the CGLE, then so is e im A{x + rj2, t + r) for any element (r)i, 772) 
of the two-torus T 2 and r £ K. Thus, we focus our study on time-periodic solutions 
of the CGLE relative to the T 2 -symmetry, namely, solutions of the CGLE that satisfy 

A(x,t) =e iv A(x + S,t + T) (1.3) 

for some (ip, S) € T 2 and T > 0. Note that after a period of time T, such a so- 
lution returns not to itself, but rather to an element in its T 2 -orbit. Such relative 
periodic solutions represent invariant three-tori in the CGLE flow. Previous studies 
of time-periodic solutions of the CGLE i|l.lfl have centered on two types of solu- 
tions. The first are single-frequency solutions of the form A(x, t) = B(x)e lujt (see, 
for example, ^3 EUJ GDI)) these are referred to as stationary solutions. The sec- 
ond type are generalized traveling waves, also called coherent structures, for which 
A(x,t) — p(x — vt)e 1<t '( x ~ vt ' > e luJt , where p and <j> are real-valued functions and ui is 
some frequency (see, for example, P HOUSED- These coherent structures reduce 
to single-frequency solutions by a change to a moving frame, that is, by a change of 
variables x — > x + vt, t — > t. The relative periodic solutions of interest in this study 
are different in that they exhibit more complicated temporal behavior than the single- 
frequency solutions just described. Studies of bifurcations to (stable) invariant two 
and three-tori for the CGLE are described in [55] and references therein. The relative 
periodic solutions we compute, which represent invariant three-tori, are obtained by 
working with fixed parameter values for the CGLE, as described below (and in more 
detail in § 13.41 and § , not via a bifurcation study. 

We employ a spectral- Galerkin procedure to compute relative periodic solutions 
of the CGLE. A solution to equation (|1.3|l is constructed and is then used to represent 
the relative periodic solutions. This ansatz consists of a Fourier series in both space 
and time, modified to include the T 2 element (ip, S) as an unknown. Substitution 
of the ansatz into the differential equation yields a system of nonlinear algebraic 
equations, the solutions of which give the desired relative periodic solutions. Using 
this procedure we were able to identify 77 distinct relative periodic solutions of the 
CGLE with parameters R = 16, v = —7, and p, — 5. These parameter values were 
chosen because for them the system exhibits temporal chaos, yet the value of R is 
such that a relatively small number of spatial modes can capture the dynamics of the 
system. 

The task of computing relative periodic solutions was undertaken as a first step in 
understanding the dynamics of the CGLE with the chosen parameters. The 77 relative 
periodic solutions are neither stationary nor coherent structures (as defined earlier), 
but two stationary solutions were also found. Aside from these two single-frequency 
solutions, we found no (truly) time-periodic solutions. For each relative periodic 
solution, the largest Lyapunov exponent ranges from 1.88 to 17.20, the sum of the 
positive Lyapunov exponents from 5.19 to 60.35, and the unstable dimension from 3 
to 8. In particular, the system displays unstable dimension variability, a phenomenon 
that a number of recent papers have indicated is central to the understanding of 
high-dimensional dynamical systems; see, for example, |331 154j . 
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In a general context, the study of periodic orbits has been a basic tool in the the- 
ory of dynamical systems for at least 100 years. Typically, there are infinitely many 
different periodic solutions embedded in a chaotic attractor. Using tools from periodic 
orbit theory, statistical averages that provide a description of the asymptotic behavior 
of a chaotic dynamical system can be approximated from the short-term dynamics of 
the (unstable) periodic solutions on an attractor of the chaotic system |13l HI] . Peri- 
odic solutions of chaotic dynamical systems are also used in techniques for controlling 
chaos. Given a set of unstable periodic solutions of a chaotic dynamical system, the 
chaos could be mitigated by applying small, time-dependent perturbations to stabilize 
trajectories onto chosen periodic solutions. This is the basis of the Ott-Grebogi-Yorke 
approach for controlling chaos; see, for example, [5U| . 

(True) periodic solutions of dynamical systems have been computed with a variety 
of methods. The use of Fourier series along with the numerical solution of systems of 
nonlinear algebraic equations is common (as in I^IEDI EZ|), but the use of multiple 
shooting algorithms for boundary value problems, which combine numerical integra- 
tion and the solution of systems of nonlinear algebraic equations, is also popular (see, 
for example, ^3 EE] )■ Furthermore, expansions other than Fourier series, for 
example, Taylor series expansions [2HJ and piecewise polynomials ^H], are also uti- 
lized. More recently, a variational method for finding (true) periodic solutions was 
proposed 

Often, the suitability of a numerical method for computing periodic solutions 
is illustrated by applying the method to low-dimensional dynamical systems (2-4 
equations). Higher dimensional systems have been considered in [TT] . [3l>| . and 
(with 16, 32, and 99 equations, respectively), where (true) periodic solutions of the 
Kuramoto-Sivashinsky equation were computed. The system studied in 1 1 1 1 has a low- 
dimensional chaotic attractor; those in and [3|>| have higher intrinsic dimension. 
In particular, the system considered in (^2) has a Lyapunov dimension of 8.8 |41| . 
with typical trajectories having four positive and one zero Lyapunov exponent. We 
also work with a high-dimensional system (62 real ordinary differential equations) 
when computing relative periodic solutions of the CGLE. The system we study has a 
Lyapunov dimension of approximately 11.52; typical trajectories in our system have 
five positive and three zero Lyapunov exponents. 

2. Symmetries and Differential Equations. In this section we review some 
standard definitions and results concerning differential equations with continuous sym- 
metries. For more information see, for example, |27| 021 QH1 OH • Properties of relative 
equilibria and relative periodic orbits are given in [2*31 l2*4l |3~H IHUl IBT] . 

2.1. Basic Definitions. Under suitable hypotheses, a vector field X, a differ- 
ential equation 

z = X(z) (2.1) 
on a space B, and an R-action or flow <f> : M x B — ► B are connected by 

d<k{z) 



x(<M*)) 



dt ' 



where <fit{z) is commonly used for cj>(t, z). The trajectory of a flow with initial condition 
zq is {<j>t(za) : t 6 R} and is also called a flow orbit or a solution. A curve j(t) is said 
to be a solution of (|2.1|l if j(t) — (f>t(zo) for some zo € B. 

Continuous symmetries are expressed by the smooth left action of a Lie group 
G on B, where we write g ■ z for the action of the element g € G on z S B. For 



4 



LOPEZ, BOYLAND, HEATH, MOSER 



simplicity of exposition we assume that G is compact and connected and acts linearly 
on the linear space B. The differential equation (|2.1[1 is said to be equivariant under 
the group action if g ■ X(z) — X(g • z) for all g € G and z € B. If <j> t is the flow of X, 
then G-equivariance of X is equivalent to 

g ■ (j>t{z) = <j> t {g ■ z) (2.2) 

for all g, z, and t, and one also speaks of the flow being G-equi variant. Since the vector 
field X is assumed to be time- independent, the equation (|2.1|l is invariant under the 
product action of G x K given by (g, r) • z(t) = g ■ z(t + r). 

The group orbit of z £ B under the G- action is G ■ z = {g ■ z : g € G}. If the 
group acting is clear from the context, this orbit will also be denoted by O z . The 
isotropy subgroup or stabilizer of z is G 2 = {g G G : g • z = z}, and the group 
orbit G • z is diffeomorphic to the quotient G/G z . As a consequence of (|2.2(l . the 
stabilizer is constant on a trajectory, that is, G^ t ( z ) = G z for all t, so one may 
speak unambiguously of the G-stabilizer of a flow trajectory or of the topological 
type of its G-orbit. In addition, if H C G is a subgroup, then its fixed subspace 
Fix(iJ) = {z E B : h ■ z = z for all h S H} is flow invariant. 

2.2. Relative Equilibria and Relative Periodic Orbits. There are two 
equivalent points of view from which one can discuss the influence of symmetries 
on the dynamics of 1)2.1(1 . The first is to consider the family of trajectories that 
arise by applying symmetry operations to a given trajectory. The second is to note 
that l|2.2[) implies that the flow on B descends to a flow on the quotient space B/G, 
and then study the dynamics of this reduced flow. 

To pursue the first point of view, note that l|2.2|l implies that given the trajectory 
with initial condition zq, the trajectory with initial condition g ■ zq can be obtained 
by applying the symmetry g to the trajectory of zq. Thus, the equivariance of the 
flow implies that trajectories come in families related by the symmetries. For a given 
zq, the union of all the members of this family is equal to the group orbit of zq under 
the joint action of T :— (G x R) on B given by (g, r) • z — <fi T (g ■ z) — g ■ <fi T (z). Thus, 
each initial condition generates a flow invariant set T ■ zq, and we may consider these 
T-orbits as the principal dynamical objects associated with the differential equation 
with symmetry 

The simplest case of a T-orbit occurs when z$ projects to a rest point of the 
reduced flow. This occurs exactly when the G-orbit of zq is flow invariant, and so 
the T-orbit of zq is equal to its G-orbit. In this case the set T ■ zo is called a relative 
equilibrium, and the point zq is said to be an element of the relative equilibrium. 
When zo is an element of a relative equilibrium, its flow trajectory is always equal to 
the action on z$ of a one-parameter subgroup o~(t) of G, that is, 

(Pt(z ) = a(t) ■ z . (2.3) 

The case of primary interest here are T-orbits that project to periodic orbits of 
the reduced flow. The sets are called relative periodic orbits. Thus, a point z$ is an 
element of a relative periodic orbit if and only if 

zq = go ■ 4>t(z q ) (2.4) 

for some go € G and T > 0. Since a relative equilibrium in the form (|2.3(l will satisfy 
(12 .4|) for go = a(t) and T — t, for any t, we explicitly include in the definition the 
requirement that a relative periodic orbit is not a relative equilibrium. The group 
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element go is called the drift, and T is the period. If Zg is an element of a relative 
periodic orbit, then its T-orbit is diffeomorphic to (G ■ z ) x S 1 . We will refer to the 
flow orbit of a point zq satisfying (|2.4|) as a (go, T) -relative periodic orbit or solution, 
and when it is clear from the context (in a slight abuse of terminology) , simply as a 
relative periodic orbit or solution. 

Note that if (|2.4|) holds, then 4>t( z o) — So ' 4>t+T(zo) for all t. Also note that 
(|2.4I) does not uniquely determine go and T: a (go, T)-relative periodic orbit is also a 
(go , fcT)-relative periodic orbit for any integer k, as well as a (gog' , T)-relative periodic 
orbit for any g' £ G Zo . However, if zo is an element of a relative periodic orbit, then 
there is a least positive T for which <pt(zo) is a (go, T)-relative periodic orbit for some 
go £ G. Thus, we adopt the convention that when a relative periodic orbit is said to 
have drift and period (go, T) the period is always this least positive T . Even with this 
convention, the drift go is still defined only up to the stabilizer G Zo . 

Finally, we remark that there is a variety of terminology used in the literature for 
relative periodic orbits and relative equilibria. In particular, the meaning of the term 
"drift" is not completely uniform, and in Hamiltonian dynamics what is called the 
drift here is called, depending on the context, a reconstruction, geometric, or other 
variety of phase. Also note that the definitions and results of this section extend easily 
to semiflows with continuous symmetries. 

3. CGLE: Symmetries and Temporal Dynamics. The CGLE in the form 
with periodic boundary conditions generates a continuous semiflow on S := 
L 2 (S' 1 ) as well as on other spaces requiring greater regularity [SJ [21 CM EH] ■ The 
semiflow on S is known to have a finite dimensional attracting set or universal attrac- 
tor as well as a Lipschitz inertial manifold |51 ll2ll22"ll26U53| . Further, on the universal 
attractor the CGLE generates a flow (i.e., solutions that are bounded and defined for 
all time), and elements of the attractor are in fact C"'' -functions |2T] . 

3.1. Symmetries and Associated ODEs. The CGLE has a number of well 
known symmetries that are central to its behavior and our study (see, for example, 
E3E]). If A(x, t) is a solution of l|l.l|l . then so are 



where rj\, 772, and r are any real numbers. Since the r\i can be treated mod 2ir, (|3.1[) 
and ll.'i.2H say that the CGLE semiflow is equivariant under the action of the two- 
torus group T 2 on S with the action on a £ S given by (771, 772) • a(x) = e lr,1 a(x + 772)- 
The symmetry (|3.3H corresponds to translation along the semiflow. The Z2-action 
expressed by (I3.4|l will be important in describing properties of specific solutions in 
§ 15.21 but in most of this paper the continuous symmetries will be our primary concern. 

Since the chosen boundary conditions for the CGLE 1)1. 1|) are periodic in x, we 
may use spatial Fourier series to transform into C :— ^ 2 (C), the space of all square- 
summable bi-infinite, complex sequences. Thus, we write 



e im A(x,t), 
A(x + r\ 2 ,t), 
A(x,t + r), 
A(-x,t), 



(3.1) 
(3.2) 
(3.3) 
(3.4) 




(3.5) 
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substitute into i|l.lfl . and obtain an infinite system of ODEs, 
da„ 



= Ra m - m 2 (l + w)a m - (1 + i/x) ^ a mi a m2 a* m!i , (3.6) 

for the complex- valued functions a m (t). 

Under this transformation the symmetries of become symmetries of (|3.6|) . 
Thus, if (a m (t)) is a solution of i|3.6[) . then so are 

(e^a m (t)), (3.7) 

(e im ^a m (t)), (3.8) 

(a m (i + r)), (3.9) 

(a_ m (i)). (3.10) 

In particular, (|3.7|) and 1)3. 8|) say that the ODEs (|3.6I) are invariant under an action 
of T 2 on C, with the action on a = (a m ) € C given by 

(Vi,m) ■ (a m ) = (c^e im ^a m ). (3.11) 

In all our computations with the CGLE we work with the spectral Galerkin 
projection obtained by fixing an even number N x , truncating the expansion (|3.5|l 
to include the indices m with —N x /2 + 1 < m < N x /2 — 1, and considering only 
the corresponding ODEs from l|3.6[) . Much accumulated theory and computation 
demonstrates that for sufficiently large N x the behavior of this truncation captures 
the essential features of the dynamics of l|l.lfl 2fl 129). 

The system of N x — 1 complex ODEs resulting from the Galerkin projection has 
the same symmetries as the infinite system. The T 2 -action is linear on the finite- 
dimensional space C^" 1 , with a group element (571 , 772) € T 2 acting via multiplication 
by the unitary matrix diag(e 1?n e lmm ) , where it is understood that the indexing is by 
increasing value of m, that is, m = —N x /2 + 1, . . . , N x /2 — 1. Henceforth this ordering 
is assumed whenever applicable. For convenience, we denote the truncated system of 
ODEs as 

^=X(a). (3.12) 

3.2. Stabilizers, Group Orbits, and Invariant Subspaces. For our discus- 
sion of the temporal dynamics and relative periodic solutions of the CGLE it will be 
useful to have a catalog of the various stabilizers, group orbits and fixed subspaces. 
These facts are well known and elementary. It is simplest to work with the T 2 -action 
on C given by (|3.11(l and then transfer the conclusions to S. 

The only point in C fixed by all of T 2 is a = 0, so henceforth we only consider 
nonzero elements and proper subgroups of T 2 . The nontrivial isotropy subgroups 
are of two classes. The first class consists of finite cyclic subgroups of T 2 gener- 
ated by elements of the form (27rfc/g, 2ir/q), for integers q and k with < q and 
< k < q. Note that k and q may have common factors. These groups will be 
denoted C(k,q). The fixed subspaces corresponding to these cyclic subgroups are 
Fix(C(/c, q)) = {(a m ) : a m = if m ^ — k mod g}, and the stabilizer of (a m ) is 
C(k,q) exactly when (a m ) & Fix(C(k,q)) and the least m > with a m ^ is 
m = q — k. The points with this stabilizer have group orbits diffeomorphic to a 
two-torus and correspond to a € S which satisfy a{x) = e l27Tk / q a(x + 2n/q). 
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The second class of stabilizers consists of one-dimensional subgroups of T 2 of the 
form {(pu,u) : u G S 1 } for integers p. These groups will be denoted D(p). The 
fixed subspaces corresponding to these one-dimensional subgroups are Fix(-D(p)) = 
{(a m ) : a m = if m ^ p}, and the stabilizer of (a m ) is D(p) exactly when (a m ) G 
F\x(D(p)). The points with this stabilizer have group orbits diffeomorphic to the 
circle and correspond to a G S of the form a(x) = a p e lpx for a fixed nonzero a p G C. 

Now recall that each subspace fixed by a subgroup of T 2 is invariant under the 
CGLE semiflow. Using the ODEs (|3.(i[) . the semiflow restricted to each Fix(D(p)) is 
easily understood. If p 2 > R, then all trajectories in Fix(D(p)) converge to zero in 
forward time. If p 2 < R, then all trajectories in Fix(D(p)) are attracted to what is 
usually called the plane wave, 

A(x,t) = ce iuJt e ipx , (3.13) 

for R = p 2 + c 2 and u> = — vp 2 — [ic 2 . 

The dynamics on the subspaces Fix(C(fc,<7)) are in general quite complicated. 
But it is worth noting that the transformation A(x,t) = q 2 B(qx,q 2 t) is a bijection 
between solutions A(x,t) of with isotropy subgroup C(0, q) and solutions B of 
(11.11) with the same v and /i but with the R coefficient equal to R/q 2 . Thus, all the 
relative periodic solutions for R — 16 described in §[5] are present, after a change of 
coordinates, in invariant subspaces for R = 16n 2 , for all integers n. 

Finally, since the focus here is on relative periodic solutions with respect to the 
T 2 -symmetry, we mention only briefly the Z2-symmetry acting as in ()3.4J) or (|3.1Q(I . 
The new class of fixed subspaces consist of elements of S that satisfy 

a(— x + 772) = a(x) or — a(— x + 772) = ct(x), (3-14) 

that is, they are even or odd with respect to reflection about 772/2. 

3.3. Relative Equilibria and Relative Periodic Solutions. One-parameter 
subgroups of T 2 can be written as (u>it,u>2t), so the relative equilibria of the CGLE 
(|1.1|) contain solutions of the form 

A(x,t) = e iuJlt a(x + uj 2 t) 

for a G S. When loi = 0, these solutions have a single temporal frequency and are 
called steady solutions. When uj\ — 0, the solution is a traveling wave, and when both 
group elements are nonzero the solution is what is called a generalized traveling wave 
or coherent structure. These have been extensively studied in the literature (see, for 
example, Q] |H EE] ) ■ 

A relative periodic orbit of the CGLE with drift (ip, S) and period T contains 
solutions that satisfy A(x, t) = e lip A(x + S,t + T) for all t. The corresponding solution 
of the system of ODEs (|3.6|l thus satisfies 

a m (t) = e iv e imS a m (t + T) (3.15) 

for all m and t. Note that the only recurrent trajectories with stabilizers D(p) are 
the plane waves (|3.13l) . which arc relative equilibria. Since by definition a relative 
periodic orbit is not a relative equilibrium, any relative periodic orbit of the CGLE 
has a T 2 -stabilizer that is either trivial or finite, and thus its group orbit is a two- 
torus. Thus, a relative periodic orbit of the CGLE is diffeomorphic to a three-torus, 
so each solution described in §E]yields an invariant three-torus in the CGLE semiflow. 
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By virtue of i|3.4|) . if A(x, t) is a solution of the CGLE, then so is A{— x, t). This 
implies that a relative periodic solution typically yields two distinct invariant three- 
tori in the CGLE semiflow. However, if the solution possesses the symmetry (|3.14[) . 
these tori are identified. 

We finally remark that the task of finding a relative equilibrium for the CGLE can 
be reduced to solving a three-dimensional system of ordinary differential equations 
with parameters (see, for example, £Q). No such simplification appears to be possible 
for finding relative periodic solutions of the CGLE. 

3.4. Dynamics for Selected Parameter Values. The parameter values R = 
16, v = — 7, and \jl = 5 were selected for our study. For these parameter values the 
numerical computation of the largest Lyapunov exponent Ai of 50 typical trajectories 
evolved from time t = to t = 600 gave a mean value of Ai ~ 5.36, with a standard 
deviation of 0.11. We obtained a value of approximately 11.52 for the Lyapunov 
(or Kaplan- Yorke) dimension, with a standard deviation of 0.04. As a consequence 
of the T 2 -symmetry, each trajectory has two Lyapunov exponents equal to zero [2], 
in addition to the zero exponent in the flow direction. These three exponents were 
included in calculation of the Lyapunov dimension. A typical trajectory has five 
positive Lyapunov exponents, and thus an unstable manifold with real dimension 5. 
The sum of the positive Lyapunov exponents (which will henceforth be called the total 
instability) was typically equal to approximately 15.71, with a standard deviation of 
0.40. 

The fairly small deviations from the mean for these dynamical quantities justify 
referring to them as "typical values" , but we use this term without any claim of the 
existence of any particular class of ergodic invariant measure. Similarly, we will speak 
of the typical trajectories as being contained in a main attractor without claiming 
that it is indecomposable in any sense. If we accept the Kaplan- Yorke and related 
conjectures, then the computed Lyapunov dimension indicates that the main attractor 
has fractal dimension equal to 11.52. The terminology "main attractor" refers to the 
closure of the typical trajectories and is distinguished from the "universal attractor" 
whose existence is proved in [5] El E21 ■ The latter, which would be called a global 
attracting set in topological dynamics, refers to the intersection of the forward images 
of a forward invariant set that is globally attracting. In particular, the universal 
attractor can contain trajectories that are not typical in the sense just described. For 
example, the Lyapunov spectra of the rest point and the plane waves, as reported 
in Appendix A, Table A. 2, are very different from those of typical trajectories and 
from those of the relative periodic solutions reported in Table A.l. 

When comparing the results here with those in the literature it is important to 
note that in addition to (|1.1|) there are other forms of the CGLE that appear in the 
literature. The most common alternative is to fix the coefficient of A as 1, 

f>A Pi 1 A 

— = A + (l + w)^-(l + ^)A\A\ 2 , (3.16) 

and have periodic boundary conditions, A(x + L,t) = A{x, t). The transformation 

A(x,t) = v r RA(v r Rx,Rt) (3.17) 

takes a solution A of (f 1 . If) to a solution A of (|3.16J) with period L = 2-k\J~R in x. In 
particular, the parameter value of R = 16 used here corresponds to solving (|3.16|) with 
L = 8tt « 25.13 and the commonly used value of L — 512 in (|3.16|) corresponds to 
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R w 6640 in The range of periods 0.02 < T < 0.46 described in §0 corresponds 

to periods 0.32 < T < 7.36 for l|3.16|l . Note also that the rescaling of time in l|3.17|l 
rescales Lyapunov exponents and so will alter the total instability of a trajectory but 
not its Lyapunov dimension or its unstable dimension. 

The value R = 16 used in this paper for the CGLE in the form is comparable 
to that of other studies of the CGLE temporal dynamics (for example, |32l I46j ) . but is 
much less than the corresponding R value used in most studies of the spatiotemporal 
behavior of the CGLE. The parameters as chosen have the advantage that the tempo- 
ral dynamics are chaotic and indeed, have a much higher intrinsic dimension than is 
usual in dynamical studies, but R is sufficiently small that a relatively small number 
of spatial Fourier modes suffices to capture the dynamics. As has been noted (for 
example in |U]), the dynamics of the CGLE at small R (or L) is quite different from 
that at larger values. As for the commonly used characterizations in terms of phase 
and defect turbulence, we note that for a typical trajectory of the CGLE at the chosen 
parameter values the modulus |^4(x,i)| frequently vanishes (defects) and the average 
phase gradient (i.e., the winding number about the origin of the image of x <— > A(x, i) 
in the complex plane as shown, for example, in FigureEH) varies somewhat irregularly 
from —2 to 2. The chosen values of v = — 7 and fi = 5 are in the Benjamin- Feir un- 
stable region 1 + fiu < 0, and the plane waves are unstable to sideband perturbations 
and have unstable dimension of 8 (see Appendix A, Table A. 2). 

4. Computing Relative Periodic Solutions. We now describe the procedure 
used to compute relative periodic solutions of the system of ODEs l|3.12|l and, thus, 
of the CGLE. Although the discussion is specifically for the problem of computing 
relative periodic solutions for the CGLE, the procedure is applicable to the general 
problem of computing relative periodic solutions of evolution equations with contin- 
uous symmetries. 

4.1. From ODEs to Nonlinear Algebraic Equations. In seeking relative 
periodic solutions a{t) of the system of ODEs (|3.12(1 . we use the ansatz 

a m (t) = e-^e-™^ a m , n e m ^ 1 (4.1) 

n 

and work with the Galerkin projection obtained by fixing an even number Nt, so 
that the summation index in l|4.1|l runs over the range —N t /2 + 1 < n < Nt/2 — 1. 
The a m (t) given by 14.1fl solve the equations l|3.15|l and hence are an appropriate 
representation for a relative periodic solution of the system of ODEs. Note that a 
(<?, T)-relative periodic solution of the system of ODEs is in fact being represented as 

a(t) = c- tL ^ T b{t/T), (4.2) 

where g = (tp, S), L g is the matrix diag(i(^ + imS), and b is periodic with period one. 
Since the action of T 2 is linear on the space C^ 31-1 , substituting (|4.2|l into the system 
of ODEs H3.12f> and using the T 2 -equivariance of X yields a system of equations for b 
and (g,T), 

from which one computes relative periodic solutions. We point out that l|4.2|l can be 
viewed as the passage into a moving frame in which the relative periodic solution is in 
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fact periodic. The moving frame is constructed using the appropriate one-parameter 
subgroup as in jSUJ. 

Indeed, in the general case of finding relative periodic orbits of a G-equivariant 
ODE, z = X(z), one adopts the ansatz z{t) = e~^/ T ■ (3(t/T), where £ is in the Lie 
algebra of G and /3(t + 1) = f3(t). Substitution into the ODE and some manipulation 
yields that z(t) is an (e*, T) -relative periodic orbit if it satisfies 

^f^?- Y e°^ = Xo ^ ( 4 - 3 ) 



T \dt 

where Y^ is the infinitesimal generator corresponding to £. The terms in l|4.3|) repre- 
sent vector fields along the closed curve (3 with Y^ o (3 corresponding to the velocity 
vector induced by passage into the moving frame. 

To be concrete, substituting l|4.1|l into the truncated system of ODEs i|3.6|) results 
in a system of nonlinear algebraic equations, 

- — - to— a m: „ = Ra m ,n - to (1 + w)a m , n 



T T T / 

, (4-4) 

mi+m2- ifi3- m \n\ J r n2—n^—n / 

for the complex Fourier coefficients {a mi „}, the drift (ip,S), and the time period T. 
Therefore, l|4.4|l is an underdetermined system of (N x — l)(N t — 1) complex equations 
in (N x — l)(Nt — 1) complex unknowns plus three real unknowns or, after splitting 
the equations into their real and imaginary parts, 2(^2 — l)(^Vt — 1) real equations 
in 2{N X — l)(JVj — 1) + 3 real unknowns. Solutions to this system of equations will 
give the desired relative periodic solutions of the truncated system of ODEs via the 
expansion (|4.1|) . 

It is convenient to separate (|4.4|l into its linear and nonlinear parts and write it 

as 

F(a,<p,S,T) = F L (a,^,5,r) + F NL (a) = 0, (4.5) 

where a is a vector with components given by the coefficients {a m .n}, Fl(o., (p, S,T) 
is a vector whose components are given by 

. / 2-nn ip S\ 2 \ 

1 1 — - j, - m j; I - if + m (1 + iv) I a m , n , 

and FNL(fl) is a vector with components given by 




Note that the components of Fnl(o) are just the coefficients in the truncated Fourier 
series expansion (in both space and time) of the function (1 + i/j,)A|^4| 2 . We remark 
that in defining the vector a (and similarly for Fl and Fnl) we are implicitly assigning 
an ordering on the coefficients {a m n \ that uniquely determines an index for the 
components of a. Henceforth, such a convention should be understood whenever 
applicable. In addition, we will use the notation in (|4.5J) to denote both the system 
of complex equations Q4.4[l and the system of equations l|4.4|l split into its real and 
imaginary parts, as it should be clear from the context which case applies. 
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4.2. Symmetries of F = 0. The symmetries ()3.7[) - (|3.10|) of the system of ODEs 
(|3.6|) induce symmetries of the system of nonlinear algebraic equations l|4.4|) . Note 
that if ({a m . n }, ip, S, T) is a solution of F = 0, then for any (ryi, 772, t) G T 3 

({^a m , n },<p,S,T), (4.6) 

({e im *a m ,„},^,S,T), (4.7) 

({e inT a m , n },^,S,T), (4.8) 

({a_ m ,„},^-S,T), (4.9) 

are also solutions. From the continuous symmetries l|4.fi|l - l|4.8|l . it follows that the set 
of solutions to F = splits into orbits O,^ v st) °^ tne s y mme try group T 3 , 

°(a, v ,s,T) : = {(Vi,m,T) ■ (a,<p, S,T) : (771, rj 2 , r) e T 3 } , 

where the action of T 3 on a point (a, S, T) is defined by 

( m , V2 ,T)-(a,<p,S,T) = ({e^e im ^e inT a m , n },<p,S,T). (4.10) 

That is, T 3 acts on a via multiplication by the matrix diag(e 1 ' 7l e lm ' 72 e mT ), and it acts 
trivially on (ip, S, T). Since the solutions of F = come in continuous families, when 
computing solutions to F = numerically we will augment the system with additional 
equations (and hence solve a system with equal numbers of equations and unknowns), 
in effect fixing a particular solution from each orbit 0^ s TV 

4.3. Kernel of Jacobian Matrix at Solutions of F = 0. Recall that, split 
into its real and imaginary parts, the system F = is an underdctcrmincd system 
of 2(N X ~ l)(N t - 1) real equations in 2(N X - l)(N t - 1) + 3 real unknowns. Thus, 
the kernel, ker(J), of the Jacobian matrix J of F is at least three-dimensional. From 
the T 3 -equivariance of F (implied by the symmetries (|4.6|I - H4.8|I ), one obtains three 
linearly independent vectors in ker( J) at a solution (a, ip, S, T) of F = 0. Such vectors 
result from a basis for the space of infinitesimal generators of the action (|4.1U|) of T 3 
on the point (a, ip, S, T), which can be obtained as 

= ({ia m , n }, 0,0,0), (4.11) 

1)1=0 

= ({ima TO , n }, 0,0,0), (4.12) 

1)2=0 

= {{in a m ,„}, 0,0,0). (4.13) 

r=0 

Split into real and imaginary parts, the above vectors V\, V2, and V3 are in ker(J) at 
a solution (a, ip, S, T) of F = 0. 

As previously pointed out, we will augment the system F = with additional 
equations and look for solutions of the augmented system. The additional equations 
should be defined so that the vectors in (|4.11l) - (|4.13() are not in the kernel of the 
augmented Jacobian matrix at a solution, since a singularity of the Jacobian is usually 
problematic for solvers of systems of nonlinear algebraic equations. One must then 
confront the question of whether, at a solution of F = 0, the vectors in l|4.11|l H4.13JI 
span ker(J), so that the Jacobian of the augmented system will be nonsingular at 
solutions. Typically this will be true, as explained next. 



v\ := 



Arji 



{{e^a m , n },<p,S,T) 



v 2 



({e im ^a m ^ n },<p,S,T) 



«3 := ({e inT a m , n },^S,T) 
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For a (<7 , T)-relative periodic solution j(t) C of a G-equivariant system of 
ODEs, where g £ G acts on C N via multiplication by a matrix N g , define the relative 
monodromy matrix M to be 

M = N go M, 

where M is the time T monodromy matrix along the solution. (For our problem 
N go = e Ls o, where L go is the matrix diag(i<^ + imS).) Thus, M is the derivative 
at 7(0) of the transformation g ■ (f>x that sends 7(0) to itself by flowing along the 
relative periodic solution for time T and then moving as dictated by the drift g . Say 
that 7 is a regular relative periodic solution if the only eigenvalues of M equal to 
1 are those coming from the continuous symmetries (including time translation) of 
the system of ODEs. Let F = denote the system of nonlinear algebraic equations 
whose solutions give (g , T)-relative periodic solutions of the system of ODEs. If 7 
is a (<7 , T)-relative periodic solution, then one can show that each eigenvalue A of 
the relative monodromy M corresponds to an eigenvalue — log(A)/T of the Jacobian 
matrix J of F at the solution to F = 6 . Now, J will have many eigenvalues that 
are not related to eigenvalues of M , but every zero eigenvalue of J does correspond 
to a 1 in the spectrum of M. Thus, if 7 is a regular relative periodic solution, then 
the kernel of the Jacobian matrix of F is spanned by the infinitesimal generators of 
the group action. Furthermore, among ODEs with a given continuous symmetry, it 
is C fc -generic for all k > 1 that all relative periodic orbits are regular |23| . 

4.4. Additional Equations. When defining additional equations to augment 
the system F = 0, we want to accomplish two objectives: first, to guarantee that 
solutions are not lost as a result of adding new equations; second, to derive conditions 
necessary so that the vectors given in l[4.11jl - l|4.13[l are not in the kernel of the Jacobian 
matrix of the augmented system at a solution. There are a variety of options available 
to accomplish these two goals. The construction we found to work best for our problem 
is to fix the sum of the arguments of some given sets A; = {a 1 1 }f =1 , 1 — 1,2, 3, of 
K nonzero coefficients by imposing the constraints 

K 

Varg(fi mi „0 = Q, (4.14) 

z ' J ' J 

3=1 

where C\, c 2 , and C3 denote some given constants. Note that in the indexing notation 
used above, the superscript I is used to indicate membership of a given coefficient 
a m i n i into one of the sets Ai, A2, or A3. To verify that imposing the constraints 14.1411 
does not result in a loss of solutions to F = 0, recall that, because of the symmetries 
(|4.6(l - (|4.8|l . if {{b m , n }, V, S, T) is a solution to F = 0, then so is {{a m< n}, <p, S, T), 
where 

n — pim pimin Jut i 

and the element (771 ,t]2,t) € T 3 is arbitrary. Therefore, the constraints 14.14f> are well 
defined if there exist some r\\,r}i,T such that (for I = 1,2,3) 

K 

^2 arg(exp(i?7i) exp(im^ 2 ) exp(m^r)6 m i i7l ^) = q 
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or, equivalently, if a solution to the linear system 



K 
K 
K 



E; 

E, 



E,™| 



E, 

















m 






T 





s) 



(4.15) 



C l - Ej ar S( fo r, 
C 2 - Ej ar S(^r, 

C 3 - Ej ar g(^ j 

exists. But one can always assign coefficients to each set A; so that the matrix in the 
above system of linear equations is nonsingular and therefore guarantee the existence 
of a unique solution to l|4.15[) . Thus, there always exists a symmetry transformation 
by which it is possible to transform any solution of F = so that it satisfies the 
constraints l|4.14(l . Hence, the constraints are well defined. 

Now let (a,(p, S,T) be a solution of F = 0. Let J denote the Jacobian matrix 
of F augmented with the equations (|4.14|) and evaluated at (a,ip,S,T), and let the 
vectors v\, $2, and V3 be as in l|4.11|) - (|4.13|) . Then it follows that 



Jvi 





K 
K 
K 



Jv 2 





E. m : 
E 7 - m 





Ej n) 
Ej n ] 

E^ 3 



where € R2(^-i)(JV t -i)_ 
least one of the sums E 7 - TO '-- 



Thus, if the sets Ai, A 2 , and A3 are nonempty, and if at 
^ and Ej n j 
will not be in the null space of J. 



n°„- are nonzero, then the vectors V\, #2, and £3 



4.5. Identifying Minimum Time Period and Relative Equilibria. It can 

happen that a computed solution (a, ip, S, T) to F = represents a relative periodic 
solution of the CGLE whose period T is not minimal. In such a case, there exists an 
integer p > 1 and rational numbers j and I so that 



f = T/p, S = (S + 2irj)/p, (p = fo> + 27rf)/p 



satisfy 



A(x + S,t + f) = A(x,t). 



(4.16) 



27r(Ti + i+jm) 



Using the expansions l|3.5|l and (|4.1|l one has that 

e^A(x + S,t + f) = ^(e-'^e-^^arn^e^e 1 

m \ n 

so ()4.16|l holds only if a m ,n — whenever p does not divide n + I + jm. The nonzero 
coefficients a m . n define a set {a m . n >}, where 

Q J m,(n+l+jm)/p = &m,n 

if p divides n+l+jm, giving a point (a, S, T) that corresponds to a relative periodic 
solution of the CGLE with minimal time period T. 

Now, a relative equilibrium A(x,t) = e luJlt a(x + L02I) of the CGLE satisfies the 
equation (|1.3H defining relative periodic solutions with (tp,S,T) = (—u>iu,—uj2U,u), 
for any u > 0. Using the expansion a(x) = E c m s lmx and such values of ((/?, S, T) in 
the ansatz (14.11) yields that for all m, a m o = c m , and when rt 7^ 0, a m>n = 0. More 
generally, a solution to F = representing a relative equilibrium will have a m , n ' 7^ 
for one index n' (not necessarily n' = 0) and a m ji = when n =/= n'. Thus, from 
the computed values of (a, <p, S,T) it is easy to detect when a solution is a relative 
equilibrium and not a relative periodic solution. 
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4.6. Starting Values for Solving F = 0. As with almost any system of nonlin- 
ear algebraic equations, the augmented system F = must be solved using a method 
that involves iteration. This requires us to supply an initial guess as a starting point 
for the iterative procedure. In general, it is very important to supply initial guesses 
that in some sense are "close" to actual solutions of the system being solved. For 
our problem, we adapted the idea of searching for close returns (or recurrences) of a 
chaotic dynamical system and sought relative close returns of the truncated system 
of ODEs to generate the required initial guesses. 

When computing solutions u(t) of a chaotic dynamical system that are truly 
periodic in time (i.e., u(t) — u(t + T) for some time period T) it is common to search 
for close returns of the system, that is, integrate the system until u(tx) « u(t 2 ) 
for some times t\ < t 2 , to within some tolerance, and use the close return as an 
initial approximation to a time-periodic solution. Since for a (g, T)-relative periodic 
solution of the truncated system of ODEs l|3.6fl one has that |a m (i)| = \a m (t + T)| 
(cf. equation (|3.15fl ). we sought sets of coefficients {a m (t)} for which 

\a m {t\)\ ~ \a m {t 2 )\ 

for some times t\ < t 2 . Given a measure S of desired closeness, the truncated system 
of ODEs H3.6[) was integrated starting with some random set {a m (0)} of coefficients 
until times t\ < t 2 were found for which 



max 



1 |fflm(*2)| 



|a m (*i)| 



: S. 



The initial value To for the time period was then set to To = ti — 1\ and a least squares 
solution over tp and S to the system 

a m {h) = e^e imS a m (t 2 ) 

was taken as an initial value go = (<po,So) for the drift. Once such a relative close 
return was identified, initial values for the Fourier coefficients {a m ,n} were generated. 
We set (cf. equation (|4.1|) ) 

e i( Vo /T )t e im(So/T )t am ^ _ ^ & m n e in ( 2 ^/ T o) t 

n 

and used the Fast Fourier Transform (FFT) algorithm to compute the initial guess 
for the coefficients {a m n }. 

4.7. Numerical Solution of F = 0. The augmented system F = was solved 
using the nonlinear least squares solver lmder from the MINPACK software pack- 
age 07] • This solver is an implementation of the Levenberg-Marquardt method 38 
1421 14*5) for solving nonlinear least squares problems. The solver lmder was chosen 
because it performed better on the CGLE problem addressed in this study than the 
other alternatives considered, namely Newton's method with a line search and 
a modification of Powell's hybrid method |471 151| for solving nonlinear equations. A 
brief discussion on the performance of these other solvers on our problem is given at 
the end of this section. 

To use the solver lmder we need to provide a routine that returns the value of 
F (augmented) evaluated at a given point (a, tp, S,T). Computing Fl(&, (p, S, T), as 
defined in (|4.5|) . is straightforward. As for Fnl(o), note that since it is a vector 
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with components given by the coefficients in the truncated (space-time) Fourier series 
expansion of (1 + hx)^4|A| 2 , one can compute it using the FFT algorithm. This is a 
standard procedure for computing nonlinear terms in spectral discretizations, see for 
example [Jy, and is briefly detailed here. Let Xj and ti denote grid-point values of 
the spatial variable x and time variable t, respectively, in physical space. To obtain 
the Fourier coefficients of A|^4| 2 , the inverse FFT algorithm was used to compute 
the values A(xj,ti), given the coefficients {a m ,n}- In order to avoid aliasing error 
when computing the convolution sums in (|4.4|l via the FFT algorithm, the coefficients 
{a m ,n\ were padded with zeros to obtain a set {a m < jn '}, where m' = — N x , . . . ,N X , 
n' = —Nt, ■ ■ ■ ,N t , and a m \n' was equal to d m ', n ' if |m'| < N x /2 and \n'\ < N t /2, 
and otherwise. (Based on the 3/2 rule for dealiasing quadratic products 8 , the 
number of modes {a m ^ n i} used are sufficient for dealiasing a cubic nonlinearity.) The 
inverse FFT was done using the coefficients {a m ' n '}. Then, the pointwise products 
A(xj, ti)\A(xj, ti)\ 2 were computed in physical space. Finally, the FFT algorithm was 
used to compute the Fourier coefficients of A|A| 2 . At this point only those coefficients 
whose index (to', n') satisfied \m'\ < N x /2 and \n'\ < N t /2 were retained. The FFTW 
software package |25) was used in our computations. 

Some points are worth mentioning regarding the additional equations (|4.14|) . 
First, pairs of indices {(mi, n'-)}jL l5 I = 1,2,3, were associated with the sets A 1; 
A2, and A3 so that the linear system in (|4.15|) was guaranteed to have a solution. 
Second, the values of c±, c%, and C3 in (|4.14|) were set to the sum of the arguments 
of the coefficients corresponding to the current approximation to a solution. That 
is, each time F augmented was evaluated, the equations 14. 14ft were satisfied exactly. 
This is appropriate since the justification for using these additional equations does 
not depend on the specific values of ci, C2, and C3. Finally, which elements were in the 
sets Ai , A 2 , and A 3 depended not only on the pairs of indices associated with them, 
but also on the value of the coefficients corresponding to the current approximation 
to a solution. Coefficients whose absolute value was below a given tolerance, usually 
10~ 5 , were excluded from the sets. 

The MINPACK software package makes use of the Jacobian matrix of the system 
being solved (or some approximation to it) and provides the option of computing 
the Jacobian using finite differences or (analytically) via a user supplied routine. We 
chose the latter option. Note that it is straightforward to compute the derivatives 
of Fl(<i, <p, S, T) with respect to ip, S, T, and the real and imaginary parts of the 
coefficients {a TOjn }. On the other hand, computing the derivatives of Fnl (a) with 
respect to the real and imaginary parts of the coefficients {a mn } directly from the 
convolution sums in (|4.4|) is potentially an error-prone process. However, one can use 
the directional derivative A 2 V* + 2\A\ 2 V of the function A| A| 2 in the direction of a 
perturbation V{x,t) to A(x,t) to compute the product of a given vector v (whose 
components are the real and imaginary parts of the Fourier coefficients {v m ,n} of 
V{x, t)) and the Jacobian matrix of Fnl evaluated at a point a. By virtue of the dis- 
cretization used, this matrix-vector product has components given by the coefficients 
in the truncated Fourier series expansion (in both space and time) of the function 
(1 + \n)(A 2 V* + 2\A\ 2 V). Hence, the matrix-vector product can be computed using a 
procedure analogous to that for computing Fnl(o), which was outlined earlier in this 
section. The difference here is that one has two sets of Fourier coefficients, {a mj „} 
and {v m ,n}j and that the pointwise products to be computed in physical space are for 
the function A 2 V* + 2\A\ 2 V. The Jacobian matrix of Fnl can thus be computed by 
substituting for v the standard basis vectors in R 2 ^- 1 )^* -1 ) . 
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As for the measure of closeness 6 used to identify relative close returns, setting 
5 > 0.5 worked well, but choosing smaller values of S usually proved unsuccessful in 
identifying relative close returns. We used the value d = 0.5, chosen experimentally, 
and the initial guesses generated from the resulting relative close returns were fre- 
quently good starting points for the solver lmder. When searching for relative close 
returns, we also specified a range [pi,P2] for the initial time period To. If 5 was 
chosen to be too large, then most relative close returns would end up having period 
T =p\. With 5 = 0.5, the search for relative close returns was always successful, but 
also resulted in relative close returns with different values within the interval [7*1,^2] 
for the initial time period To. We remark that the use of random initial data (i.e., 
not resulting from a relative close return), even with the spectral profile of a typical 
trajectory of the truncated system of ODEs, was unsuccessful in generating solutions 
with any of the nonlinear equations solvers. 

Initially, the system of nonlinear algebraic equations was solved using values of 
N x = 32 and N t = 48 (giving a total of 2,917 real variables). For this value of N x , 
the decay in the spatial spectra of the relative close returns was around five orders 
of magnitude. This proved to be sufficient in order to have well defined solutions in 
the sense that once a solution was found, continuing from it with N x > 32 posed no 
challenge for the nonlinear equations solver and in a very small number of iterations 
it converged to the same solution. The value Nt — 48 was chosen to allow for a 
relatively large range of time periods. If the solver lmder converged to a solution 
of F = 0, then we verified whether the solution had minimal time period using the 
criteria from § 14.51 If not, a new set of Fourier coefficients was defined and the values 
of T, S, and if modified in order to obtain a solution with minimal time period. Then 
the values of N x and Nt were increased, if necessary, to have a good decay (around 
eight or more orders of magnitude) in the spectra of the solution. At this point, using 
Newton's method to solve the system of nonlinear algebraic equations worked well: 
only 2-4 iterations were required for convergence. Finally, the solution to F = was 
used as an initial condition and the truncated system of ODEs (|3.6(l was integrated 
to verify that we indeed had a relative periodic solution of the system of ODEs (and 
thus of the CGLE). 

As for the other solvers considered, the performance of Newton's method with a 
line search |16| was, in general, unsatisfactory for our problem. The typical behavior 
was that, as the iterations progressed, the line search parameter became very small and 
no significant change in the iterates nor reduction in the value of ||F||| was achieved. 
The MINPACK solver for systems of nonlinear algebraic equations, hybr j , was also 
tried on our problem. It is an implementation of a modification of Powell's hybrid 
method |47l 15 1| for solving nonlinear equations. The solver hybr j often terminated 
without converging to a solution of F = nor to a minimum of F T F for which F / 0. 

5. The Relative Periodic Solutions. 

5.1. Dynamical Properties. Using the procedure described in §0J 77 distinct 
relative periodic solutions were found for the CGLE (|l.l(l with the parameter values 
R = 16, v = —7, and \x = 5. The solutions display a wide variety of dynamical 
and spatial properties, some of which are summarized in Appendix A, Table A.l, and 
illustrated in Figures fo~71 and 15751 In Table A.l the solutions are sorted by increasing 
period, with the periods in the range 0.02 < T < 0.46. The solutions are given an 
identifying number in the left column corresponding to their position in the list. 

The Lyapunov exponents of the solutions were computed from the eigenvalues of 
the relative monodromy matrix defined in § 14.31 As a consequence of the continuous 
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symmetries and the flow direction, each relative periodic solution has three zero expo- 
nents PJ. This was confirmed numerically. Figure luHT a) shows the largest Lyapunov 
exponent of each solution plotted against the corresponding time period, while Fig- 
ure EHb) shows the total instability. The horizontal dotted line on each graph gives 
the value of the corresponding dynamical quantity computed for a typical trajectory, 
which is assumed to be characteristic of the main attractor (cf. § 13.4(1 . 
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Fig. 5.1. (a) Largest Lyapunov exponent vs. period and (b) total instability vs. period for 
relative periodic solutions. Mean value for typical solutions denoted by horizontal dotted line. 

Figure 15.21 shows the unstable dimension of each solution plotted against the 
period. As noted in § 13.31 each relative periodic solution represents an invariant 
three-torus in the CGLE semiflow, and so, for example, a relative periodic solution 
with five positive Lyapunov exponents corresponds to a flow invariant three-torus 
with a 5-dimcnsional unstable manifold. The unstable dimension variability of the 
system is clearly illustrated in Figure 15.21 We have been unable to associate the 
various values of unstable dimension with specific regions of the main attractor, but 
the high dimensionality makes a clear geometric picture of the attractor very difficult 
to obtain. 
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Fig. 5.2. Unstable dimension vs. period for relative periodic solutions. Value for typical solu- 
tions denoted by horizontal dotted line. 

As seen in Table A.l, the solutions possess a wide variety of space shifts S and 
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phase shifts ip. A relative periodic solution can be the "k th root" of true periodic 
solution if A(x, kT) — A(x, 0) for some integer k. This happens exactly when kip 
and kS are both integer multiples of 2n and corresponds to the dynamical situation 
of a rational flow (every orbit is periodic) on the invariant three-torus. To test for 
this possibility we computed the best rational approximation to each (p/(2n) and 
S/(2ir) using the MATLAB function rationalize. For no solutions did the rational 
approximations of (p/(2ir) and S/ {2n) share the same denominator, and in most cases 
the denominators of the rational approximates were on the order of several hundred. 
These results indicate that, at least for moderate size k, the computed relative periodic 
solutions are not k th roots of true periodic solutions and that the dynamics on the 
invariant three-tori are either an irrational flow or else a periodic flow with long period. 

The semilog plots of the temporal and spatial power spectra of all the computed 
relative periodic solutions are approximately linear, indicating near exponential decay 
in these spectra. Figure I5T41 in § 15.21 shows plots of the spectra for solution 14. The 
exponential decay reflects the fact that solutions on the CGLE attractor are in 
time [52) and space [21) . In addition, the temporal spectra indicate that none of the 
77 computed relative periodic solutions are coherent structures. 

5.2. Properties of Specific Solutions. To give an idea of the character of the 
solutions we describe a few in some detail. The first relative periodic solution found, 
identified as solution 14 in Table A.l, is the least unstable among our collection of 
relative periodic solutions. The time evolution of solution 14, represented as curves 
on the plane with coordinates defined by the real and imaginary parts of A(x, t) at 
different times in the interval [0,T], is displayed in Figure [^31 The nontrivial time 
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Fig. 5.3. Imaginary vs. real part of solution 14 at times t = 0, T/3, 2T/3, T . Solution has phase 
shift tp = 2.61. Click here for movie of time evolution. 

behavior of the solution (i.e., the excitation of multiple temporal frequency modes) 
is apparent from these curves. For single-frequency solutions A(x,t) = B(x)e lu)t and 
generalized traveling waves A(x,t) = p(x — ti^e 1 ^ - "*^ 1 "*, where to is some single 
frequency, plots of this kind at different times would show, except for a rotation, the 
same curve. It is clear, then, that solution 14 is not of either of these types. Note also 
that the curve at time t = T differs from that at time t = because of the rotation of 
the complex field A(x, 0) — > e lv> A(x + S,T). The spectra of solution 14 are displayed 
in Figure I5~4l where one can also observe the multiple modes active in the temporal 
spectrum of the solution. 

Another of the solutions with total instability less than typical, listed as solution 
16 in Table A.l, has period and drift close to those of solution 14. However, solution 
16 has stabilizer (7(1,2) (cf. §E2J, that is, it satisfies A(x,t) = -A(x + n,t). This 
property is apparent from the symmetry in each of the curves representing the time 
evolution of solution 16, displayed in Figure 15.51 We will come back to solutions 
14 and 16 in § 15.31 where an apparent relationship between these and several other 
solutions is explored. 
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Fig. 5.4. Spectra of solution 14- 




Fig. 5.5. Imaginary vs. real part of solution 16 at times t = 0, T/3, 2T/3, T . Solution has phase 
shift if: = 2.72. It satisfies A(x,t) = — A(x + 7T,t). Click here for movie of time evolution. 

Contour plots of the absolute value of a relative periodic solution provide a way of 
visualizing periodicity in time coupled with the shift S in space. Figure IST^l shows one 
such plot for one of the most unstable among the relative periodic solutions found, 
listed as solution 20 in Table A.l. The absolute value of the solution has been plotted 
over two space and four time intervals. Moving vertically (in time) and horizontally 
(in space) allows us to observe periodicity in time, after taking into account the shift 
S in space. From the contours of \A\ one can also locate points in space-time where 
the magnitude of A vanishes (the darkest shaded regions (in blue) in Figure^))). The 
time evolution of solution 20 is shown in Figure 15.71 To view the time evolution of 
several other relative periodic solutions, click on the following red links: solution 8, 
solution 24, and solution 26. 
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Fig. 5.6. Contours of \A(x,t)\ for solution 20, x £ [— 27r,2ir], t £ [0, 4T]. Solution has space 
shift S = 1.27. 
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Fig. 5.7. Imaginary vs. real part of solution 20 at times f = 0, T/3, 2T/3, T . Solution has phase 
shift tp = 5.81. Click here for movie of time evolution. 

Finally, several other relative periodic solutions have nontrivial T 2 -stabilizers and 
some are even or odd about certain points in the spatial domain. For instance, 
Figure 15.81 displays contour plots of the real and imaginary parts and the absolute 
value of a relative periodic solution having all these properties. The solution, identified 
as solution 15 in Table A.l, has total instability close to that of a typical trajectory. 
The drift and period of solution 15 are close to those of solutions 14 and 16, but the 
unstable dimension of each of these solutions is different. Like solution 16, solution 
15 has stabilizer C(l, 2). In addition, solution 15 is odd about x — it and even about 
x = 7r/2,37r/2. Since solution 15 has stabilizer C(l,2), it is also a relative periodic 
solution with period T and drift ((p, S) = (tp + ir, S + ir). The space shift S of solution 
15 is (numerically) equal to w, so it follows that its absolute value is (truly) periodic 
in time, with time period T. All these properties of solution 15 can be observed from 
the contour plots in Figure l5~8l 
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Fig. 5.8. Contour plots for solution 15, x € [0,2tt], t € [0, 2T]. Solution satisfies A(x,t) = 
— A(x + ir,t) and \A(x,t)\ = \A(x,t + T)|. It is odd about n for x £ [0,27r], even about ir/2 for 
x £ [0,7r], and even about 3iz/2 for x £ [7r,27r]. 
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The other relative periodic solutions with nontrivial T 2 -stabilizers are solutions 3 
and 9, which have stabilizers C(l,2) and C(l,3), respectively. Solutions 4 and 10 are 
even about x — tt and have a space shift S (numerically) equal to tt. The latter results 
in the absolute value of each of these solutions being periodic in time with period 2T. 
Such properties of solutions 3, 4, 9, and 10 can be seen from Figures f5.9H5.lll We 
remark that solutions 3, 9, and 10 are more unstable than typical trajectories, whereas 
solution 4 is among the least unstable of the solutions found. 
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Fig. 5.9. Imaginary vs. real part of solution 3 at times t = 0, T/2, T . Solution satisfies 
A(x, t) = — A(x + 7T, t). Click here for movie of time evolution. 




Re(A(«,0)) Re(A(K,T,'21) R,.(A..Ti) 



Fig. 5.10. Imaginary vs. real part of solution 9 at times t = 0,T/2,T. Solution satisfies 
A(x,t) = e i27T / 3 A(x + 2n/3,i). Click here for movie of time evolution. 




5.3. Instability and a Family of Solutions. There are several features of 
the instability of the collection of relative periodic solutions that are apparent from 
Figures 15.11 and 15.21 the largest Lyapunov exponent and total instability of most 
of the relative periodic solutions are larger than those of a typical trajectory, and 
in contrast, the unstable dimensions are nearly evenly distributed about the typical 
value. In addition, the periods of the relative periodic solutions with instability less 
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than typical in Figure lSTl are clustered around what appear to be multiples of a single 
fundamental period. 

More specifically, if Tf = 0.054, then the seven relative periodic solutions whose 
largest Lyapunov exponent and total instability are less than typical have periods 
close to one of Tf, 2Tf, 4Tf, and 6!/. As will be discussed below, the drifts (<p, S) 
of the corresponding solutions are also approximate integer multiples. This suggests 
that these least unstable relative periodic solutions arose via a sequence of period- 
multiplying bifurcations and that the resulting collection of solutions are members of a 
low instability (and thus high traffic) region of the main attractor, which is structured 
around a three-torus that projects in the reduced flow to a loop with the fundamen- 
tal period Tf. Further evidence for this proposed structure is the observed repeated 
appearance in typical trajectories of sequences of patterns from the low instability so- 
lutions. The role of R in the equations (|3.6|) makes it clear that, in general, instability 
will increase with increasing R. Thus, low instability orbits could be continuations 
of solutions present for smaller values of R. This conjectured attractor structure and 
bifurcation scheme should certainly be investigated in further detail. 

In the initial phases of this investigation we have found a collection of solutions 
that clearly fit together in a dynamical family. The family contains four of the seven 
solutions with total instability less than typical and, in fact, consists of four of the 
five solutions with least total instability. The family splits into the pair 14/16, which 
have periods approximately 2Tf, and the pair 58/59, which have periods approxi- 
mately 67/. The elements of a pair have very similar (ip, S,T), and all these values 
for the pair 58/59 are close to three times those of pair 14/16. Solutions 14 and 16 are 
distinguished from each other by their unstable dimensions of 4 and 5, respectively. 
Solutions 58 and 59 both have unstable dimension of 3, but they differ in their un- 
stable Floquet spectrum (i.e., the eigenvalues of the relative monodromy matrix with 
magnitude greater than one). In this portion of the spectrum, both solutions have a 
conjugate pair of complex eigenvalues, but the remaining real eigenvalue is positive for 
solution 58 and negative for solution 57. Thus, they have opposite Lefschetz indices 
as fixed points of their Poincare map and hence the available information is consistent 
with solutions 58 and 59 arising together in an equivariant period-tripling bifurcation. 

All four of the solutions in the family have similar appearing pattern evolution, 
and indeed, this is one of the features that brought them to our attention. In partic- 
ular, solution 16 has stabilizer C(l, 2) and thus satisfies A(x + ir,t) — — A(x, t), while 
the other three solutions approximately have this symmetry. To investigate these 
"near symmetries", recall from § 13.21 that a solution has a (7(1,2) stabilizer exactly 
when its spatial Fourier coefficients (a m ) satisfy a m — for all even m. Thus, for 
each of the solutions we plotted the magnitude of the even and odd parts of the power 
spectrum using 



The results are shown in Figure I5.12f a^ , which clearly shows the strong similarities 
in the evolution of d\(t) and the differences in do(t) among the members of the fam- 
ily. Furthermore, the magnitude of di(t) is significantly larger than that of do(t), 
confirming the presence of the approximate stabilizer. As further confirmation of the 
similarities in the evolution of the coefficients <Z2fc+i(£) for members of the family, 
we show in Figure l5.12f b) the projection of the solutions onto the three-dimensional 
space with coordinates (|a_x(t)|, |ai(t)|, |a3(i)|) (these are the three a,2k+i(t) with 
largest magnitude). We confirmed comparable similarities among the evolution of 
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(a) (b) 

Fig. 5.12. (a) Magnitude of even (do(t)) and odd (d\ (I) ) parts of power spectrum for solutions 
14, 16, 58, 59. (b) Projection of time evolution onto three-dimensional space with coordinates 
(\a-i{t)\,\ ai (t)\,\a 3 {t)\). 

the other |a2/c+i(^)|- Figure I5\l2f b) also shows the pair 58/59 traversing the other 
solutions three times. 

To put these findings in a more geometric context, we consider a direct sum 
decomposition of the space C into its odd and even pieces. More precisely, let Fo = 
{(a m ) : a m = if to is even} and F\ = {(a m ) : a m = if m is odd}. Thus, 
F = Fix(C*(l,2)), Fi = Fix(C(0,2)), and C = F ffi F\ . Furt her, f or i = 0,1, the 
L 2 -distance of (a m (t)) to Fj is di(t), and so Figures 15. 12f al and 15.121^ ) demonstrate 
that the projections of the family of solutions onto F are all approximately equal, 
and that the time evolution of the distances to Fi are all approximately the same. 
In other words, the family of solutions comes very close to lying in the generalized 
cylinder A x Fi, where A c Fo is the full relative periodic orbit of solution 16 (i.e., 
its invariant three-torus), and so the bifurcations that (presumably) created these 
solutions respected this cylinder. Since the values of d\ are consistently larger than 
those of do, the solutions are evolving near this cylinder reasonably close to the Fi 
subspace. 

Under the Fourier transform, the decomposition C = F © Fi corresponds to a 
decomposition S = Fq © F[, where for a S S the projections onto Fq and F[ are given 
by (l/2)(a(x) — a(x + n)) and (l/2)(a(a;) + a(x + ir)), respectively. Thus, within the 
family the evolutions of (1/2) (A(x, t) — A(x + ir, t)) are all very similar while those of 
(l/2)(A(x,t) +A(x + <K,t)) differ. 

The existence in the family of a pair with T « 2Tf and a pair with T w 6Tf 
raises the question of solutions with T « Tf and T « ATf. The question of solutions 
with T w Tf asks whether the pair 14/16 is itself the result of a period doubling 
bifurcation. Solution 4 has T ~ Tf and it, along with the family, constitutes the five 
solutions with the least total instability. Using its drift and following it for two of 
its time periods, solution 4 satisfies ^4(2;, 0) = exp(5.868i)At(:r, 0.108). On the other 
hand, shifting the drift of solution 16 by its stabilizer yields for one time period the 
very similar 

A 16 (x, 0) = exp(5.857i)^i 6 (x, 0.114). (5.1) 
Thus, in terms of its drift, period, and total instability solution 4 fits in with the 



24 



LOPEZ, BOYLAND, HEATH, MOSER 



family, but it lacks the near symmetries that characterize the family and is located 
nowhere near the family's cylinder. The situation for solutions with T w 4Tf is some- 
what similar. After reducing mod 2tt, the drift and period of solutions 30 and 31 are 
approximately twice those of solution 16 as given in (|5.1|) . Thus, while the seven solu- 
tions with total instability less than typical have drifts and periods in (approximate) 
integer multiples, only four of the solutions show clear signs of being related by well 
understood bifurcations. 

5.4. Weighted Averages of Observables. As mentioned in the introduc- 
tion, the unstable periodic orbits embedded in a chaotic attractor can be used to 
approximate statistical averages that characterize the long-term dynamics of the sys- 
tem ^H^J[1}2]. In an initial effort to investigate the usefulness of the relative periodic 
solutions in this regard, we use a weighted averaging scheme proposed in |62| for sys- 
tems for which a symbolic dynamics is not known. 

A given quantity $ t yp attached to a typical trajectory is approximated by 

N 

$ typ «< $ > N := ^p^i, (5.2) 

»=i 

where 

l 

_ SL, 

$i is the value of the quantity $ computed for the z-th periodic solution, and SL^ 
is the sum of the positive Lyapunov exponents of the i-th periodic solution. The 
value 1/SLj approximates the fraction of time that a typical trajectory spends in a 
neighborhood of the i-th periodic solution. 

The quantities $ we compute are: the largest Lyapunov exponent Ai; the Lya- 
punov dimension Dl = j + (Xh=i where Ai > A2 > . . . > A m are 
the Lyapunov exponents and j is the largest integer such that — ®> the 
time average of the Ginzburg-Landau energy functional, E = ^ J^(—R\A\ 2 + 
\dA/dx\ 2 + \ A\ 4 /2) dxdt; the time average of the L 2 -norm of the vector field that de- 
fines the CGLE, L 2 -VF = ± / Q T '(f£* \ RA+(l+iu)d 2 A/dx 2 -(l+i^)A\A\ 2 \ 2 dx)^ 2 dt; 
and the time average of the L 2 -norm of a solution A(x, t) of the CGLE, L 2 -A = 
T So (So \ A \ 2 da;) 1/2 di. Table IO summarizes the results obtained using the N = 77 
relative periodic solutions found. (We remark that the rest point and the plane 
waves (cf. Section I3.4|l were not included in the computation of the weighted aver- 
ages.) From the relative errors (fourth column in Table I5TT|) one could say that, except 
for the value of Ai, the approximations obtained are quite good. 
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While these initial results are promising, a detailed investigation of the predictive 
value of the relative periodic solutions is warranted. For example, other averaging 
schemes which do not require knowledge of symbolic dynamics, such as those pre- 
sented in [171 118} . should also be considered. Another question concerns the role that 
unstable dimension variability played in the computation of the trajectory averages 
<&typ since, as discussed in j^i], convergence to erroneous values for statistical averages 
computed from long-time simulations is a possibility in the presence of unstable dimen- 
sion variability. As for the relative periodic solutions, since they are computed using 
the expansions ((3.51) and (|4.1|) (as opposed to integrating the system of ODEs (I3.fc>|) ^ 
we do not expect unstable dimension variability to have had an adverse effect in the 
computation of the values $i . 

5.5. The Iterative Solver and Dynamical Properties of Solutions. With 
the rather large range of stability properties exhibited by our collection of relative 
periodic solutions, it is interesting to consider whether these properties impact the 
ability of the nonlinear equations solver to converge to solutions of F = 0. The 77 
relative periodic solutions listed in Table A.l resulted from 350 runs performed using 
the nonlinear least squares solver lmder and relative close returns (with a measure 
of closeness of 5 = 0.5) as starting values. In total, 146 of the 350 runs resulted in 
convergence to solutions of F = (this count includes convergence to the 77 distinct 
relative periodic solutions, runs that converged again to these solutions, and two 
runs that yielded single-frequency solutions of the CGLE), and 204 runs converged 
to minima of F T F for which F / 0. For those runs that converged to solutions of 
F = 0, the final value of ||F||2 was on the order of 10 -9 or smaller. The total number 
of iterations taken by the solver lmder ranged between 10 and 300. 

For solutions with periods in the range 0.40 < T < 0.46, the value N t = 64 
(instead of Nt = 48; see § 14. 7f) was required to obtain well defined solutions, in the 
sense that the decay in the temporal spectra of the solutions found was around five 
orders of magnitude. For a generic compact attractor there will be finitely many 
relative periodic orbits in any finite range of periods. Of the total number of runs 
performed with the solver lmder, the last 55 yielded three new solutions. This gives 
around a 5.4% return on the number of new solutions found from the last set of 
runs, which suggests that we were approaching the point where no new solutions 
with periods less than 0.46 could be found with our choice of solution method. This 
assertion is difficult to prove, however, based only on empirical testing. 

Because we solve the equations (|4.4f) using starting values derived from relative 
close returns of the system (I3.6f) coming from random initial data, the method we use is 
biased towards finding relative periodic solutions on the main attractor. Nonetheless, 
there is no reason to expect the iterative solver to respect the dynamics of the CGLE, 
so solutions disjoint from the main attractor are certainly possible. However, there do 
seem to be quantities associated with the CGLE dynamics that are, in certain cases, 
connected to the solver behavior. 

Solution 14 (shown in Figure 15.3(1 was by far the most frequently found by the 
solver, with convergence to it 25 times out of the 146 runs that converged to solutions. 
It was also the first solution found as well as being the least unstable solution by both 
measures shown in Figure 15.11 The instability of a relative periodic solution can be 
connected (cf. § 14.3(1 with the eigenvalues of the Jacobian matrix of F at the root, 
which in turn is, at least hcuristically, connected with the eigenvalues of the Hessian 
of the Newton method iterator at the root. Thus, low instability could be connected 
with the presence of a large basin of attraction for the iterative solver. On the other 
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hand, again heuristically, under the dynamics of the CGLE the least unstable relative 
periodic solution would be the one that is easiest for other orbits to approach, and 
thus low instability would tend to indicate that the solution is in a high traffic area of 
the attractor, that is, in an area with a large concentration of invariant measure. As 
noted in § 15.31 observations of the time evolution of typical solutions do reveal their 
occasional return to a sequence of patterns similar to those manifested by solution 14. 
This means that the starting values we used for the solver are more likely to be in a 
neighborhood of solution 14, which would favor convergence to that solution. It is not 
clear which mechanism, a large basin or an abundance of nearby starting values, is 
most responsible for the frequency of occurrence of solution 14 among the computed 
relative periodic solutions. 

While solution 14 was the most favored solution, there are other instances of 
multiple convergence. However, these do not, in general, have the same correlation 
with the total instability. On one hand, the second most frequently found solution was 
solution 16, with convergence to it 9 times. This solution does have low instability and 
is in fact paired with solution 14 in the family described in § 15.31 On the other hand, 
the next two most frequently found solutions, solutions 3 and 11, had 6 convergences 
each, but they do not have low instability. In addition, solution 58 has total instability 
only slightly more than that of solution 14, but the solver found it just once. So while 
we can, in certain cases, connect the solver behavior to the CGLE dynamics, it is, 
in general, extremely difficult to understand the behavior of iterative solvers in high 
dimensions. 

6. Discussion and Conclusions. We have described a method of finding rela- 
tive periodic solutions for differential equations with continuous symmetries and have 
demonstrated its utility by finding numerous relative periodic solutions for the CGLE, 
a nonlinear evolution equation with widespread physical and mathematical interest. 
The various purposes to which periodic orbits are usually put, such as dynamical 
averaging and control, can be adapted to using relative periodic orbits. In our inves- 
tigation we found an abundance of relative periodic solutions and an apparent lack of 
true periodic solutions. This indicates that in the presence of continuous symmetries 
the use of relative periodic orbits is an attractive option. In addition, the computa- 
tional burden of finding relative periodic orbits is not significantly higher than that of 
finding true periodic orbits. It would be of interest to adapt other methods of finding 
periodic orbits to find relative periodic orbits and systematically compare the various 
methods. 

The most computationally challenging aspect of the numerical procedure is the 
solution of a system of nonlinear equations for the drift, period, and Fourier coefficients 
of a relative periodic solution. The use of starting values derived from relative close 
returns and the choice of solver were critical to our success in identifying solutions. 
The size of the system considered here, between 2,917 and 3,909 real variables, is 
at the upper limit of feasible computations at workstation scale. In particular, the 
solution of linear systems with a (dense) Jacobian as coefficient matrix, which was 
required as part of the procedure for solving the nonlinear system, could be done 
using direct methods. For systems with higher dimension than that considered here, 
one must eventually incorporate the iterative solution of linear systems, along with 
preconditioning, and the use of parallel computations. 

As noted in the introduction, the computation of the relative periodic solutions 
was undertaken as a first step in understanding of the dynamics of the CGLE with the 
chosen parameters. This prompts two questions, one specific and the other more gen- 
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eral. First, is the computed set of relative periodic solutions sufficient to capture the 
total dynamics in any reasonable sense, and second, is the computation of a collection 
of relative periodic solutions a good starting point for analyzing and understanding 
the dynamics of a moderately high-dimensional dynamical system with continuous 
symmetries? 

The general method of analyzing dynamics using a collection of periodic solutions 
has been investigated in many ways, for many years. A given single periodic orbit 
will usually not be typical with respect to any globally supported ergodic measure, 
but under hypotheses that guarantee the existence of a Bowen-Margulis measure (see 
§ 20.1 in HO]), the total collection of periodic orbits yields invariant measures that 
are equidistributed, and thus the collection of periodic orbits captures the statistics of 
the system. While these strong hypotheses can never be confirmed in practice, there 
is a fair amount of empirical evidence that a sufficiently large collection of periodic 
orbits can capture the global dynamics, but it is usually impossible a priori to decide 
what is "sufficiently large" . 

The next issue is what part of the dynamics one can hope to capture with a 
collection of periodic orbits. A typical moderate dimensional chaotic system will 
display a great deal of dynamical variety, and there is no reason to expect a single, 
indecomposable attractor. As noted in § 15.51 the method we use is biased towards 
finding solutions in the closure of the set of typical trajectories. Thus, the appropriate 
question about the collection of relative periodic solutions is in what sense it predicts 
the behavior of a typical trajectory. 

We have begun this investigation, and the preliminary results are promising. As 
reported in § 15.41 typical values of the Lyapunov dimension and various functionals 
were accurately computed using the 77 relative periodic solutions and the weighted 
averaging scheme proposed in |fi2| . We have had no success in finding symbolic dynam- 
ics, and preliminary attempts to use the sequence of patterns given by the various 
relative periodic solutions to build up the pattern evolution of a typical trajectory 
have been unsuccessful. Thus, the question of the predictive value of our collection of 
relative periodic solutions remains open. 

The second question raised above we can answer with an unqualified "yes" . Dy- 
namical systems with even moderate dynamical dimension present significant chal- 
lenges not adequately dealt with using the usual low dimensional theory. The family 
of relative periodic solutions has allowed us to begin unlocking the intricacies of a 
flow with Lyapunov dimension of 11.52 (including the flow and symmetry directions). 
The detailed study of the dynamics of systems with comparable intrinsic dimension 
is still a rare occurrence in the literature. In this regard we remark that it is com- 
mon in certain segments of the literature to use the terms "high-dimensional chaos" 
or "hyperchaotic" for any system with unstable dimension larger than one, and the 
unstable dimension in our system is typically 5 and varies from 3 to 8 among the 
relative periodic solutions. 

We have found the relative periodic solutions invaluable as a starting point of 
dynamical analysis for a variety of reasons. First, one can be quite sure of the ac- 
curacy of the numerical computations because the solutions can be verified both by 
integrating the differential equations for a short time interval, [0, T], and by con- 
firming that the drift, period, and Fourier coefficients are solutions of the system of 
nonlinear algebraic equations. Next, each computed relative periodic solution can be 
studied in isolation as one piece of the dynamics, as one part of the pattern evolution. 
Finally, the observations of the individual pieces leads to hypotheses and conjectures 
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about the whole picture. These can be tested using the collection of relative periodic 
solutions as a sample space and then the results compared to typical trajectories. 

As an illustration of this process we briefly list some of the lines of inquiry that 
have been opened by our study of the collection of relative periodic solutions. Is there 
a correlation between spatial aspects of the patterns of a solution and its temporal 
dynamics, for example, between winding numbers and total instability? How does the 
unstable dimension variability of the collection of relative periodic solutions contribute 
to the existence of zero finite-time Lyapunov exponents for the system and other 
consequences of non-hyperbolicity? What are the global dynamics of our system, 
for example, how many nontrivial indecomposable attractors are there and how do 
the invariant fixed symmetry planes fit in with the attractor(s)? What role do the 
coherent structures or relative equilibria play in the overall temporal dynamics of 
the system; are they contained in the main attractor, or are they isolated and thus 
dynamically nontypical? Finally, how are the relative periodic solutions and overall 
dynamics created via bifurcations as R increases from zero to 16? 

Nonlinear dynamics has made significant progress in the study of systems with 
low intrinsic dimension, but the understanding of even moderate dimensional systems 
still presents significant challenges. We have demonstrated that finding a collection of 
relative periodic orbits is both computationally feasible and conceptually valuable, and 
is thus a valuable tool in the important task of understanding moderate dimensional 
dynamical systems with continuous symmetries. 

Appendix A: Properties of Solutions. Table A.l lists the period T and drift 
(ip, S) for each of the 77 relative periodic solutions found. The solutions are listed by 
increasing value of T and have been assigned an identifying number (first column). 
Also listed is the dimension of the unstable manifold of each solution, along with the 
largest Lyapunov exponent and total instability (which is given by the sum of the 
positive Lyapunov exponents). 

While most of the relative periodic solutions have trivial T 2 -stabilizers, there are 
a few exceptions. Solutions 3, 15, and 16 have stabilizer C(l,2) (cf. § I3.2fl . and 
therefore satisfy A(x,t) — — A(x + it, t). Solution 9 has stabilizer C(l,3); it satisfies 
A(x,t) = e l27T / 3 A(x + 27r/3, i). In addition, solutions 4, 10, and 15 are even about 
certain points in the spatial domain and solution 15 is odd about certain points in 
the spatial domain. 



Id 


T 


9 


S 


Unstable 
Dimension 


Largest 
Lyapunov 
Exponent 


Total 
Instability 


1 


0.02333 


5.36228 


3.85441 


4 


8.85592 


28.53243 


2 


0.05394 


2.88491 


3.09563 


5 


5.00321 


21.31125 


3 


0.05394 


0.00111 


3.97095 


5 


11.00430 


30.30192 


4 


0.05403 


2.93438 


3.14159 


4 


3.57030 


9.56834 


5 


0.05471 


4.60939 


1.45371 


5 


6.45679 


27.60085 


6 


0.05561 


4.51654 


4.70614 


5 


7.76507 


25.39928 


7 


0.06080 


0.24364 


2.38876 


4 


11.31755 


36.02733 


8 


0.08255 


4.79599 


3.08241 


5 


10.08510 


39.81868 


9 


0.08748 


0.28765 


2.44316 


6 


9.72294 


36.31152 


10 


0.08950 


5.02514 


3.14159 


3 


10.72851 


27.85965 


11 


0.10458 


2.60234 


3.17193 


4 


5.83393 


18.52259 


12 


0.10797 


2.65754 


3.12097 


3 


6.59523 


19.55427 


13 


0.11065 


6.05532 


0.00321 


6 


5.10800 


19.17581 



Table A.l 
Properties of relative periodic solutions. 
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Id 


T 


H> 


S 


Unstable 
Dimension 


Largest 
Lyapunov 
Exponent 


Total 
Instability 


14 


0.11282 


2.60639 


3.10577 


4 


1.88230 


5.18967 


15 


0.11461 


2.25004 


3.14159 


3 


6.90750 


15.73770 


16 


0.11492 


2.71671 


3.14159 


5 


3.20606 


10.12255 


17 


0.13560 


1.82110 


2.04598 


6 


10.56790 


42.15027 


18 


0.13730 


0.73763 


0.75086 


6 


13.54288 


53.74900 


19 


0.14037 


1.87712 


5.73712 


7 


10.57182 


40.07190 


20 


0.14916 


5.81089 


1.26773 


5 


17.19835 


50.08455 


21 


0.15160 


1.94728 


5.40241 


5 


6.75385 


19.66934 


22 


0.16098 


1.95738 


3.85354 


5 


10.61619 


44.38707 


23 


0.16263 


2.56407 


5.09938 


6 


11.43816 


35.46420 


24 


0.16971 


5.48891 


3.93522 


8 


7.86379 


30.95111 


25 


0.17500 


3.80256 


0.98842 


6 


10.36722 


44.23884 


26 


0.18249 


3.40328 


6.14235 


5 


10.19717 


36.97288 


27 


0.18317 


3.61256 


5.24560 


6 


12.29563 


42.02092 


28 


0.20146 


1.90508 


2.41894 


6 


14.46845 


59.36105 


29 


0.21038 


0.70319 


1.90349 


5 


14.85499 


60.35458 


30 


0.21178 


5.16441 


0.07994 


4 


5.87695 


14.88747 


31 


0.21766 


5.61399 


0.00002 


5 


3.70755 


13.08357 


32 


0.22045 


5.15235 


0.06622 


3 


8.28842 


16.37539 


33 


0.22482 


4.41470 


1.10856 


7 


13.46647 


47.71035 


34 


0.24430 


5.53858 


3.60540 


6 


11.40020 


41.36390 


35 


0.24598 


0.88868 


4.47924 


5 


8.78079 


35.50550 


36 


0.24761 


2.24375 


2.32864 


6 


15.96544 


59.60573 


37 


0.25234 


4.06727 


2.73577 


5 


7.27440 


28.49855 


38 


0.25374 


1.47983 


2.26736 


6 


10.89903 


43.88382 


39 


0.26014 


5.12178 


2.60985 


4 


5.39730 


16.59190 


40 


0.26623 


3.07369 


3.27604 


6 


11.19883 


50.85904 


41 


0.27286 


3.69060 


3.55692 


6 


9.88356 


37.26567 


42 


0.27504 


2.57465 


3.21456 


5 


5.92789 


18.84754 


43 


0.27679 


2.90953 


3.10889 


4 


6.77859 


24.98180 


44 


0.28744 


4.84490 


1.76031 


5 


7.34615 


23.42338 


45 


0.28803 


2.74448 


1.37428 


5 


8.98143 


31.75536 


46 


0.29746 


4.15104 


4.09627 


6 


6.98488 


22.53727 


47 


0.29751 


4.32786 


1.58049 


5 


6.18589 


20.33417 


48 


0.30045 


4.39783 


4.62959 


4 


5.77264 


18.52750 


49 


0.30392 


3.72216 


5.49134 


6 


8.10462 


31.10158 


50 


0.30492 


2.74872 


3.31013 


5 


6.89564 


21.04823 


51 


0.30695 


4.31730 


5.12453 


4 


9.02700 


24.32585 


52 


0.31343 


0.72363 


1.98525 


5 


8.41168 


25.64814 


53 


0.31517 


1.87776 


1.69240 


5 


7.79446 


24.72635 


54 


0.31924 


5.68308 


2.25711 


6 


11.97732 


46.26127 


55 


0.32350 


2.21845 


1.26058 


5 


8.44337 


31.00286 


56 


0.33150 


2.42184 


3.37971 


3 


6.90072 


15.94646 


57 


0.33217 


4.50138 


1.88596 


5 


6.86741 


19.11005 


58 


0.34035 


1.72898 


3.15957 


3 


3.92346 


5.99484 


59 


0.34148 


1.74043 


3.18504 


3 


2.56575 


6.92418 


60 


0.34825 


6.00790 


5.07384 


6 


10.50662 


47.95424 


61 


0.34892 


2.45502 


3.11195 


4 


5.98084 


18.44338 


62 


0.35567 


0.76100 


0.97929 


6 


8.70424 


25.31565 


63 


0.36395 


1.11216 


5.51746 


5 


4.75266 


16.16998 


64 


0.36621 


3.81362 


2.39863 


6 


5.79290 


27.12694 


65 


0.37306 


5.25206 


1.42035 


5 


9.96569 


30.85402 


66 


0.37413 


1.14509 


1.19984 


5 


8.30154 


18.70350 


67 


0.37619 


1.78357 


0.93027 


5 


6.29217 


16.68876 


68 


0.37893 


1.44760 


0.83458 


4 


6.92225 


20.55719 



Table A.l 

Properties of relative periodic solutions, continued. 
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Id 


T 


V 


S 


Unstable 
Dimension 


Largest 
Lyapunov 
Exponent 


Total 
Instability 


69 


0.38444 


1.83435 


1.06175 


5 


7.33543 


20.40029 


70 


0.38599 


1.21873 


1.00009 


4 


7.66270 


17.44702 


71 


0.39128 


2.38456 


5.34383 


4 


5.78585 


18.67097 


72 


0.39675 


4.79784 


3.35390 


6 


9.55893 


30.02312 


73 


0.40557 


3.21458 


6.24425 


4 


6.80844 


22.92050 


74 


0.42108 


0.44223 


4.10959 


4 


10.55221 


19.44133 


75 


0.42447 


2.27338 


3.25431 


5 


7.58996 


26.84563 


76 


0.44044 


3.32603 


2.31594 


6 


7.92846 


26.18850 


77 


0.45840 


5.73253 


6.25633 


5 


8.01104 


25.60890 



Table A.l 

Properties of relative periodic solutions, continued. 



Table A. 2 lists properties of the solution A(x,t) = 0, which corresponds to ld=0, 
the plane waves A(x, t) = a P: ±ie ±luJt e lpx , which correspond to the identifying numbers 
1-4, and the two single-frequency solutions A(x, t) = B(x)e~ luJt (Id=5,6) that were 
computed numerically. Solution 6 has stabilizer C(l,2) and is even and odd about 
about certain points in the spatial domain. Solution 5 is even about certain points in 
the spatial domain. 



Id 


T 




Unstable 
Dimension 


Largest 
Lyapunov 
Exponent 


Total 
Instability 









14 


16.0000 


168.0000 


1 


0.07854 


80.0000 


8 


54.7935 


320.2662 


2 


0.09240 


68.0000 


8 


51.7819 


298.3522 


3 


0.19635 


32.0000 


8 


41.2843 


220.1846 


4 


0.22440 


28.0000 


8 


26.0693 


124.4494 


5 


0.12517 


50.1968 


10 


15.8062 


109.0126 


6 


0.12609 


49.8320 


7 


15.0947 


61.4864 



Table A. 2 

Properties of solution A{x,t) = (ld=0), plane waves A(x,t) = a P: ±ie^ zlujt e lpx (Id=l,2,3,4 corre- 
spond to p=0,l,2,3, respectively), and single-frequency solutions A(x,t) = B(x)e~ m,t (Id=5,6). 
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